_ AD-ALS4 552 DEVELOPMENT OF REDUCED ORDER MODELS FOR CONTROL SYSTEM a/72.. 
DESIGN USING THE Sa ai PROGRAM(U) oe POSTGRADUATE 
SCHOOL MONTEREY CA W NELSON DEC 8 
UNCLASSIFIED F/G 972 


———* 


{ 
{ 
| 
4 


FEE 


| eg 


FREE 
EEE 


FEREEE 


— 
° 
— 
—E 

& 

fe 


—— 
S— 
— 
—— 


S 
== 
a 
= 
o fil 


MICROCOPY RESOLUTION TEST CHART 
NATIONAL BUREAU OF STANDARDS-1963-A 


! 


NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


AD-A154 552 


THESIS 


DEVELOPMENT OF REDUCED ORDER MODELS 
FOR CONTROL SYSTEM DESIGN USING 
THE OPTSYSX PROGRAM 
by 
Stanley William Nelson 


December 1984 


DTIC FILE COPY 


Thesis Advisor D.J. Collins 


CG: for public release; distribution is unlimited 


ELECTE X 
JUN5 1965 


i aa ea 


UNCLASSIFIED _ 
SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


___ REPORT DOCUMENTATION PAGE 
+ REPORT NUMBER 2. GOVT ACCESSION BO «eRECIPIENT’S CATALOG NUMBER 
a 


4 TITLE (and Subtitte) 5. TYPE OF REPORT & PERIOD COVERED 


Development of Reduced Order Models Master's Thesis 
for Control System Design Using December 1984 


7. AUTHOR(a) 6. 


Stanley William Nelson 


9. PERFORMING ORGANIZATION NAME ANO ADORESS RGJECT, TASK 


December 1984 
13. NUMBER OF PAGES 


1S. SECURITY CLASS. (of thie report) 


a 
z 
c 


Naval Postgraduate School 
Monterey, California 93943 


11, CONTROLLING OFFICE NAME ANO ADORESS 


Naval Postgraduate School 
Monterey, California 93943 


) MONITORING AGENCY NAME & ADORESS(!f different from Controlling Office) 


Unclassified 


Sa. OECL ASSIFICATION/ DOWNGRADING 
SCHEOULE 


Approved for public release; distribution is unlimited. 


16. DISTRIBUTION STATEMENT (of this Report) 


17. DISTRIBUTION STATEMENT (of the abstract entered in Block 20, If different from Report) 


. SUPPLEMENTARY NOTES 


. KEY WORDS (Continue on reverse side If necessary and identify by block number) 


Reduced Order Modelling 
Model Reduction for Control Systems 
F100 Turbofan Engine 

X-29A Aircraft Control System 


20. ABSTRACT (Continue on reverse side if necessary and identify by block number) 


The modern controls engineer is often faced with designing 
a system which is characterized by a large number of first order 
differential equations. It is highly desirable and sometimes 
necessary that such complex systems be reduced for analysis, 
synthesis and implementation into a physical system. It is the 
intent of this thesis to present a mathematical procedure and 


DD ipan 98 1473 so EOI TION OF 1 NOV 65 18 OBSOLETE 


5. N 0102-LF-014- 6601 1 SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


5 = aes = 7 a a ae, 
a DT fe a a a Me a 8 Me We ee Met ate ee Me em 


UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE (When Deta Entered) 


* a 
“6-8 @ 
oie H 


and computer software based upon this procedure which enable 
- the control engineer to formulate a reduced order model of a 


m large order system. 
As examples, two large order systems are analyzed: a fi 
- sixteenth order model of the F100 turbofan engine and a ninety- 
wO eighth order model of the X-29A aircraft control system. 
re 
: aes 
; Ce 
po 
< 
oe 
oe = 
x: is 
a *. 


= x 
a) 9 woe 
as wee 
““ cons 
aa S/N 0102- LF 014-6601 aD 
es UNC I ASS [ E I ED ; ww! 

7 2 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) ‘ 

: Lr 


Se Ty Fg Plame a, a a 
we ta teh ea 


cou Pr ere ae 
ie Se Se A A a ES 


Approved for public release; distribution is unlimited 


ar Development of Reduced Order Models 
for Control Systea Design Using 
the OPTSYSX Program 


by 


. erentey W. Nelson 
Lieutenant Commander, United States Navy 
-S., University of Kansas, 1972 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN AERCNAUTICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 


December 1984 sot 

4 ‘\ TattN 

Author: A on W. Mike a as 
SE ey W. Nelson Sr, noes 

Approved by:______ Oy (th ahs eet 
ees Pay ee Iins, Thesis Advisor poco 


ie & 
“Woanshr tis Setond header eae 


tr Fe PALIT 


To ASF. Platzér, Chairian, 
Department of Aeronautics 


Accession For 
NTIS GRA&l a panne 
DTIC TAB 

Unannounced 
Justification 


ws 
ae 
or 
RC 
3 ee” 
ot 
‘se 
ae San 
oy 
é we te et tee ioip tg tie tole ete ae te wee pt ee i. Spiele” iwi gia toa tue Bagh gi? ates Mis ate Sele wie ts Seen ate er ween Bat 
ee a ee ge ne eo ee ia te Ma ae ee woe le ag al a ae ad ea ae a Oa ae ae 


: z " yi * we. ia ve une ane Br) Leow oar Se oe ge ty, 
Ne te ie Oe ee tN at ee a at a mee aM We PM ee oe, SME Oa Se oa. : 


ABSTRACT 


» The modern controls eagineer is often faced with 
designing a system which is characterized by a large nunber 
of first order differential equations. It is highly desir- 
able and sometimes necessary that such complex systems be 
reduced for analysis, synthesis and implementation intoa 
physical control system. It is the intent of this thesis to 
present a mathematical procedure and computer scftware Lased 
upon this procedure which enakle the control engineer to 
construct reduced order models. 

As examples, two large order systems are analyzed: a 
sixteenth order model of the F100 turbofan engine and a 
ninety-eighth order model of the {X-29A aircraft control 
system. Li¢le-lin [x , OPTSYSX 


a 
! 


Level, Kakrver ) 


pat 


“2m lb. i i 7 
i faeie trl forces. Fo Fem / | ; 


atewte 


he 
° 


Tit. 


IV. 


aes hye MEE MGR fan Sg ea eee Cee Oe, ene aCe oe hy 


TABLE OF CCNTENTS 


INTRODUCTION eo e e e e oof e e e e e e @ e e e 


THE OPTSYSX COMPUTER PROGRAM .« 2. « 2 © « « «© « 
A. PROGRAM OVERVIEW . 2 2 «© © © 2 2 © © © © 
Be. OPTSYSX MODIFICATIONS . . «© « « © © © «© « 


ANALYTICAL DEVELOPMENT OF A REDUCED ORDER 
MODEL e e e e e e e es e e e e oe e e ° e e e e 


THE OPTRED PROGRAM « « « « « © © «© © © e « «© eo 
A. PROGRAM OVERVIEW . -« « « «© © © « © © © «© «© 
B. PROGRAM OPERATION . 2.2 «© « « «© « © « «@ © @ 
1. Program Composition and Flow ..... 
2. User Protection Features . ..... -e 
3. Large Order Systems .......««« 
C. INTERFACE WITH OTHER ANALYSIS PROGRAMS .. 
D. EXAMPLES OF ORDER RELUCTION ......e .« 
1. Example of a Sixteenth Order System . 

2. &Example of a Ninety-eighth Order 
System . se © 2 ew we we we we ew ew sw 


CONCLUSIONS AND RECOMMENLTATIONS . «2... « « « 
A. CONCLUSIONS . « «© « « « © «© © «© © «© «© @ 
B. RECOMMENDATIONS .. 2. « «© «© © © © © @ © « 
1. Comparison of System Time Response .. 
2. Interface with Other Analysis Programs 
3. Alternate Order Reduction Techniques . 


APPENDIX A: THE OPTSYSX PROGRAM AS MODIFIED BY THIS 
THESIS 2 e e e e s e e e e es e es s e e s 


11 


14 
14 
16 


18 


23 
23 
24 
24 
27 
28 
2? 
29 
30 


35 


49 
49 
49 
49 
50 
50 


51 


Se ee a ei ace re the ee, Pe eee See eel SLE ES ene ae ana Gee Loe eC eS Nea oe ce Ow ie cet aap 
ee ae aa a Te 


APPENDIX Bs: THE OPTRED PROGRAM ~ . « © «© © « «2 © © © «= 141 Re 


APPENDIX Cs THE OPTRED EXEC « «© e e e# @ ee e@ @ . <« e e 153 P pes 


LIST OF REFERENCES e es e e oe e e e e e eo e 4 e e es s es 154 wy 
ws 
BIBL IOGRAPHY e e e e e e es e e e e s e e es e e e es e e 155 eae 


INITIAL DISTRIBUTION LIST . . 2. 2. «© © © «© «@ «© «2 © «© © «= 156 


~ 

ws 
Varseeu 
ut 

8 
a 
® 
be 
ere 
aiens 
ete 
weer) 
. 
‘ 


® . 
eat 
“et Ah 
wos 
a 


4.1 
4.2 
4.3 
a.4 
4.5 
4.6 
4.7 
4.8 
4.9 
4.10 


LIST OF FIGURES 


Full System Bode Plot for the F100 Engine . 
Reduced System Bode Plot for the F100 Engine 


X~29A 
X-29A 
X~-29A 
X-29A 
X-29A 
X~29A 
X-29A 
X-29A 


Full System Without Notch Filter . 
Reduced System Withcut Notch Filter 
Full System Without Notch Filter . 
Reduced System Withcut Notch Filter 
Full System With Notch Filter ... 
Reduced System With Notch Filter . 
Full System With Notch Filter ... 
Reduced System With Notch Filter . 


av ats 
‘= te ne 


vu A Ww 
ft 


wO = 


ee Ma et eT ee 8 = 
ee date, See ee 


leat Se Se Sk Oe ew. 


SYMBOLS 


State (Ns,Ns) or Output (No,No) Weighting Matrix 
Control (Nc,Nc) Weighting Matrix 

Control Gain Matrix (ic,Ns) 

Control (No,Nc) or Noise (No,Ng) Feedforward 
Matrix 

Reduced Control Feedficrward Matrix (No,Nc) 
Open-Loop Dynamics Matrix (Ns,Ns) 

Reduced Open-loop Dynamics Matrix (Nr,Nr) 
Control Distribution Matrix (\Ns,Nc) 

Reduced Control Distribution Matrix (Nr,Nc) 
State Disturbance Distribution Matrix (Ns,Ng) 
Measurement Scaling Matrix (No,Ns) 

Reduced Measurement Scaling Matrix (No,Nr) 
Estimator Gain Matrix (Ns, No) 

Number of Controls 

Number of Process Noise Sources 

Number of States 

Number of Observations or Measurements 


Reduced Model Number of States 

White Process Noise Ccvariance Matrix (Ng,Ng) 
White Meas. Noise Covariance Matrix (No,No) 
Steady-State Covariance Matrix of Control (Nc,Nc) 
Control Vector (Nc,1) 

White Measurement Noise Vector (No,1), With 


ap ae Le eR tek 
we iManage Ore ee 
to. eae) aes 
24 ¢ oft 

oot, Pare hae 


. 
“i 


‘ 


zero Mean and Covariance Matrix R 


eal ghia 


Pa ee ee ree 


White Process Noise Vector (Ng,1), with 


a mers 
fake 
aft 


Zero Mean and Covariance Matrix Q 
Constant Disturbance Vector (Ng, 1) 
State Vector (Ns, 1) 

Derivative of State Vector (Ns,1) 


Dee a Pen el AST colar a mes - “at WE ELS hn Oa ae a ae ak eal 


xe 
xe 


RESET Ser Ma eee ee eS ae ee ee ea alt 


Estimate of State Vector (Ns, 1) 

Derivative of Estimate of State Vector (Ns,1) 
Output Vector (No, 1) 

Measurement Vector (Nc, 1) 


eo . 
aye * 
>. "ig 
aes . ” 
oie + * 

OS act ia 
wet 

' "a " 
aA 
. =a a”. 
~ 6 
ete” 

. whe 
% ay aI 


ms 
as . 
ae 
ae 
evs 
Cre 
re) 
veeey 
Sere 
eta he 
Seay 
origi 
io er) 
~ 


wats 


6 ES APS 28k F) cay teh ET Mg Foe ee RL RR aR ey tac eige Ads Met a UNS SRS re Rg ale ee ma a 
Pg IRS SS gay Ducat el Wied Cee gE Te eS 


ACKNOWLEDGESENTS 


I would like to express my Sincere appreciation to 
Professor D.J. Collins whose assistance ani encouragement 
enabled this research. I would also like to thank Professor 
D.G. Williams, Director of the W.R. Church Computer Center, 
for his assistance and cooperation in providing the addi- 
tional computer memory allocaticn which was essential to the 
completion of this work. 

This thesis is dedicated to my wife Gayle, my daughter 
Kristin and my son Scott. Without their constant love, 
support and understanding this work would not have been 
possible. 


10 


A er eer ee ree ee 
nr en ee ee 


en of at ae, Mae 
me oe Daee tse te 

Carat ans ek Peer Sriikty W3001, 6 wae BS 
et we at ek 


I. INTRODUCTION 


Modern control methods allow the control engineer to 
develop a control system to regulate every parameter of a 
physical systen. He has been aided by computer simulations 
of non-linear systems and linear approximations of those 
systems. The desired contrcl system can be developed 
parallel to or even ahead of the physical system which is to 
be controlied. 

In .many design applications areas, the engineer is 
confronted with designing a ccntrol system which can be 
represented analytically by a very large number of first 
order, linear, time invariant differential equations. This 
is especially true with highly complex and aero-elastic 
aircraft and space vehicles. The analytical model is a 
prime ingredient in the controller design process for any 
design technigue. 

Many practical limitations exist, such as computational 
requirements, that require the order (number of equations) 
of this complex model to be reduced for analysis, synthesis 
and actual implementation of the control system. Some method 
is required to analyze these dynamical models and establish 
Simpler systems which include elements critical to the 


desired control function. Without such simplification, 
application of design procedures can result in highly 
coap nd parameter-sensitive controlled systems. With a 
reduced 2x model, regulator synthesis procedures become 


more intui. 2 and far less sensitive to parareter varia- 
tion. The ‘educed order model is far easier to handle in 
formulating he control system but must characterize the 
physical system swith sufficient accuracy such that perforn- 
ance objectives for the controlled physical system can be 
met by designing control laws from the reduced order model. 
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Both file and interactive data entry is accomplished 


by the main program and by the primary subroutine, REDUCX. 
The actual computation of the reduced order model is 
performed by subroutine REDUCX which utilizes subroutines MI 
and MAMULT for its matrix inversion and amultiplication 
reguirenents. The remaining Subroutines are utilized 
throughout OPTRED for interactive entry of numerical -and 
character string data. The frogram reguires no external 
subroutines for its execution. 

The main program begins by presenting the user with 
a brief operational description of OPTRED. At this point, 
the user is given the option of having the program's general 
results repeated to the screen cr printed to a listing file. 
OPTRED creates a permanent data file which contains the 
reduced data regardless of the user's option here. The 
program then reads full System flags and parameters, 
excluding matrices, from the OPIMAT data file. A brief 
description of the full system under investigation is then 
presented as an indication to the user that the desired data 
is being input. This description includes the order of the 
system, the number of controls and observations and whether 
or not a system "D" matrix will be input. The latter is 
indicated either by " "or "net" in the following example. 
The program then prcmpts the user for entry of the desired 
order of the reduced systen. 

The main prcgram presentation as it appears during 
execution: 

OPTRED 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
"FT NG" "HY, AND "D™ (IF INPUT) MATRICES. 


THE FULL SYSTEM MATRICES MUST BE READ FROM A FILE NAMED 
“OPTMAT DATA" AS CREATED BY THE OPTSYSX PROGRAM. 


YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) 
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1) The compiled text file of OPTRED 
2) ##The OPTMAT data file 
3) The OPTRED EXEC file 
The first two required files are self explanatory. The 
OPTRED EXEC file sets up the required read/write file defi- 
nitions and calls the necessary library functions for execu- 
tion using the IBM 3033 System. This EXEC allows the user 
to execute the OPTRED program by simply entering the werd 
“OPTRED" at the terminal. The prog.am then continues ina 
user-interactive mode until completion as described in the 
next section. The OPTRED EXEC file is listed in Appendix C. 


B. PROGRAM OPERATION 


OPTRED was written to satisfy two specific cbjectives. 
The first is to provide a methcd to assist in the analysis 
of actual large order control systems. The second is to 
create a program to be used as an instructional tool for 
students studying controls related problems. 

To accomplish these objectives the emphasis while 
creating OPTRED was to make the program as user oriented as 
possible by eliminating data entry ambiguities and providing 
features which prevent the user from other inadvertant or 
invalid program entries. A serious effort was also made to 
Minimize the amount of virtual memory required to execute 
this frogran. 


1. Program Composition and Flow 


OPTRED is composed cf one main program and seven 
subroutines. The program can be divided into three kasic 
catagories: 

1) File Data Input 
2) Interactive Data Input 
3) Calculation 
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A. PROGRAM OVERVIEW 


OPTRED is an interactive, double precision FORTRAN 
program which uses a specifically formatted data file to 
input full system data. The required matrix inputs fron 
this file are the full system "F", ug, “He and "pf" 
matrices. The system flags and parameters which OPTRED must 
also read for program execution are Ns, Nc, No, and IFDFW. 
The OPTSYSX program creates OPTMAT data which is the source 
of this reguired data. Through interactive input, the user 
specifies all other data which is required to create reduced 
model data. 

OPTRED also reads from the CPTMAT file the noise, feed- 
back control and the output and control cost weighting 
matrices. These are not reguired inputs but only input to 
the program to facilitate creation of a reduced order data 
file which is identical in format to the OPTMAT data file. 
The final result of OPTRED is the creation of OPTMATR, a 
data file which contains the reduced system "F", "G", "q" 
and "D" matrices and is formatted for re-entry into the 
OPTSYSX frogran. 

It is possible for the user to construct a data file 
named OPTMAT to provide input data to OPTRED. However, only 
in rare applications would the full system not be analyzed 
prior to the development of the reduced nodel. OPTSYSX is 
therefore the simplest and most preferred method of data 
generation. 

To execute the OPTRED prcgranm, a Minimum of three 
separate files must be immediately accessible to the user: 
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Y= ar*x, + Dr*0 (3.24) : 
where 5 . 
Fro= Py, - Fyo*Foo-t*F 51 (3.25) : 
Gr = G, - PF), 4F,,-1*6, (3.26) 
Hr =H, - H,*F,,-*F,, (3.27) . 
Dr = D - #L*F,,-14*6, (3.28) 


This is the mathematical prccedure which OPTRED uses in 
computing reduced order models. The complete OPTRED listing 
is included as Appendix B. 
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Gy = 
Gt = eae dees (3.16) oe 
Py) a. 
HY = [H, 1 Ho] (3.17) = 
diye. x ate 
eR. 
If the states in X, truly model the full order systen, * 
then the subvector X, will be essentially zero and oe 
— * * * ‘ - 
x, Pui XK + Fi x, + Gy U (3.18) io 
= * * x » a 
0 Foy xX + F.5 xX, + G. U (3.19) aa 
Solving for X5 in Equation 3.19, one has 
= =F, .-1% * -1%*G * 
X, Foo Pou x + Fo2 G. U (3.20) 
This is incorporated into Equation 3.18 and one 
= - *Po —1k + -F. *F  -1% . ee 
x, Pia Pi Po? Pot) (G) Pio Foo G4)u (3.21) ene 
and the output eguation becomes ee 
= 
= -H_*P -1% —H_*F -1% 2 ree 
Y (H, i, Poo Posy + (TL H. Poo 6)U (2.22) as 


The reduced order model can now be written as 


X = Po*X + 3r* ; ee 
= Prek, + 3rtu (3.23) . 
a 
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X* = R#X = TZ" = {---- (3.9) 


This X* vector is the reordered state vector, where the 


elements contained in the subvector xo are those states 


associated with the designers choice of eigenvalues while 
the subvector Xo contains all other states. 
The linear system may now be rewritten 


X* = R¥F*R-1%4X" + R*G*T (3-10) 
and 
Y = H*¥R-1*X" + D*O (3.11) 


Now define 


Ft 


R*F #R-1 (3.12) 


G* = R*G (3.13) 


tee} 
ul 


H*R-1 (3.14) 


and 
' 
F ' Ff 


ll i 12 
Fi = nn anne (3.15) 
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The eigenvalues which are tc be included in the initial 
reduced order model are determined through modal analysis. 
The eigenvalue matrix can be re-ordered by constructing a 
Matrix Q, which has the value 1 in the position corre- 
sponding to the eigenvalue to ke included (column) and the 
re-ordered position (row) of that eigenvalue. For example, 
if eigenvalue 20 of the original A matrix was to be 
re-ordered to the first position, a 1 would be placed in the 
(1,20) position of the Q matrix. Now, the reordered matrix 


At = C¥A (3.4) 


is incorporated in the linear model 


Q*Z = A'*Z + Q¥T-14G40 (3.5) 
letting 
Z* = Qez (3.6) 
thus 
ZY = A'kQ-18Z4 + QeT-14G8U (3.7) 
and 
X = T*Q-1t#z0 (3.8) 


A matrix R can be constructed such that 
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III. ANALYTICAL DEVELOPMENT OF A REDUCED ORDER MODEL 

There are numerous methods of order reduction for 
control systems as discussed by Enns {Ref. 8] of Stanford 
University. This thesis uses the theory discussed by DeHoff 
and Hall [Ref. 9] which presents a method of reducing the 
order of the linear system based on dominant states chosen 
by the designer. This dominance is determined by modal 
analysis of the eigensystem under investigation. Dougherty 
[Ref. 10] utilized this same methodology in the development 
of a computational model used in the analysis of the F100 
turbofan engine. A summary of this technique is provided 
below. 

The method of reduction involves finding the eigenvector 
transform matrix, T, such that 


X*T = THA (3.1) 


where A is the diagonal eigenvalue matrix. Then, by 
defining an alternate state vector, Z, where 


X = T4Z (3.2) 


the original system may be re-written in modal coordinates 
as 


Z = A*Z + T-1%Gey (3.3) 
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ure, sgt, wat, and "GAM" matrices were complemented by the 
addition of routines which also save the "D" matrix. These 
hew routines are analagous to existing "save" routines. In 
addition to these routines a new flag, ISAD, was added to 
the existing flags ISAF, ISAG, ISAH, and IGAM. This new 
flag provides the same function in "saving" matrices as do 
the existing flags. Por program continuity, it was also 
necessary to add ISAD to subreutine INNER, which is the 
driving subroutine in the OPTSYSX progran. 

Subroutines RDMATF, RDMAT, and WRTMAT also. required 
minor changes to accomodate the "D" matrix. These changes 
consisted of the addition of "D" to the subroutine parameter 
lists and the addition of read or write statements for 
either re-entry of "D" into OPISYSX or for writing "D" to 
the OPTMAT data file. The existing flag, IFDFW, was also 
added to the parameter lists of these routines. This flag 
indicates whether or not a "D" matrix will be input by the 
user. OPTRED uses this flag in its computational methods in 
the same manner as does OPTSYSX. Appendix A lists’ the 


OPTSYSX program as modified by this thesis. 
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additions to that program's size would cause user difficulty 
ey in obtaining sufficient memory for program execution. It is 
for this reason that Diel's OPICALC, Laptas' OPTGRAPH and 
~ the program developed by this thesis have been developed as 
independent programs which depend upon various data created 
by OPTSYSX. The minor modifications this thesis introduces 
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to the OPTSYSX program are not detrinental to its original 
capabilities, 


Be OFPTSYSX MODIFICATIONS 


This thesis is ccnfronted with the problem of analyzing 
very large control system nodels. Therefore, the additions 
of Hoden's subroutine SETUP and Diel's subroutines RDMATF, 
RDMAT, and WRIMAT were crucial to system analysis. These 
routines make multiple computer analysis runs’ possible 
without painstaking re-entry of system parameters and 
matrices. In particular, the system "F", "G", "H", and 
"Gam" are some oz the matrices which are stored for re-entry 
into OPTSYSX or stored as data for use by other analysis 
programs. Subroutine WRTMAT is responsible for the storage 
of these matrices in the data file "OPTMAT". It is this 
‘4 file which OPTRED utilizes tec input all pertinent full 
system data. 

In many control applications, a control feed-forward 
matrix ("D") is present in the system model. Therefore, 
this thesis makes provisions within OPTSYSX which enable the 
user to save the "D" matrix in the same manner by which the 
other system matrices are saved. Subroutine SETUP, which 
enables the user to input system matrices from file data, 
. was changed to enable the user to input the "D" matrix. 
This consisted of adding "D" to the routine's parameter list 
and also added a read statement to input that matrix froma 
Ce data file. Present routines in OPYSYSX which "save" the 


a 


I 


a 
= 


open loop transfer function 

CHI*|s(I} - (FJ-2*(6° (2.5) 
closed loop noise transfer function 

CH}*[s(I} - (FJ[-1*{Gan] (2.6) 


compensator transfer function from measurement to input 


[C}*[s[I FL F]+(G }*(C ]+(K J*#[H]]-2*(K] (2.7) 

where 

x = state vector (Ns X 1) 

x = derivative of the state vector (Ns X 1) 

u = control vector (Nc X 1) 

y = output vector (No X 1) 

Zz = measurement vector (No X 1) 

w = white process noise vector (Ng X 1) 

[F] = open-loof dynamics or plant matrix (Ns X Ns) 

{G] = control distributicn matrix (Ns X Nc) 

(Gam} = state disturbance distribution matrix (Ns X Ng) 

(Hj = measurement distribution matrix (No X Ns) 

[Dj = control feed-forward matrix (No X Nc) 

[Cc] = control gain matrix (Nc X Ns) 

[1] = identity matrix (Ns X Ns) 

OPTSYSX is an extremely large and complex program which 
contains over 4000 lines of code. Its existing standard 
dimensions will accomodate a thirty-second order system with 
up to 10 controls and 10 observations or measurements. To 
execute even this version of OPTSYSX the user requires one 
megabyte of virtual storage. Any further significant 

15 
at eet seg ES een ee ACTER ROR eC EE ae peste ee 14 Pe Soon eee. a 


ees 


7 


wrere 
. 


rs 
f 


Py, 


s ; eine 
aa * + 
Cee a ae he a sata 
Pa eee * 
s ee ae ee eh 1 fe 
7 ola te Se te Me te as 


we 


7 Cee 
reel ee! 


tate 


a 


uf 


II. THE OPTSYSX COMPUTER PROGRAM a, On 


A. PROGRAM OVERVIEW 


Lr v * ee ad 


OFTSYSX is a double precision, interactive FORTRAN 
program which employs modern ccntrol theory analysis tech- 
nigues. It is developed to be compiled and executed by the = 
Naval Postgraduate School's IBM 3033 System 360/370. Its . 
primary capabilities include the calculation of the open 2 
loop eigensysten, and the fixed closed loop system; the L. 
synthesis of optimal regulators or filters; the power i. 
spectral density, and modal distribution computations. is 

The fundamental system equations used by the OPTSYSX e 
program for its computations are of the state variable forn. b 
The basic system equations are: - 


system model 
x = (F]*x +[(G]#u + [GAM]*w (2.1) 


output equation 


TO EE 


y = [H}*x + (D)j*u (2.2) 


measurement egquaticn 


atl A 


z= [(HJ*x + v (2.3) 


estimator equation 


xe = [F]*xe + [G]*u + [KJ*(z - [H]*xe) (2.4) , 

c 

14 : 

L 

pi pltie uate ate gree te ee ae tah ee eh Seas Te ae Ee Se te gata lee ee cathe ein tain loom ie tee tien tad tear uarsge ie Sgar te termi tpeete San oie 


systems and the corresponding reduced godels which OPTRED 


generates are presented and analyzed. Those systems include 
a sixteenth order model of the .F100 turbofan engine anda 
ninety-eighth order model of the X-29A aircraft. Complete 
program listings for the OPTSYSX program and the OPTRED 
program are included as Appendices A and B respectively. 
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The intent of this thesis is to present the mathematical 


basis for creating such a reduced order model and to develop 
the actual computer software which creates that reduced 
model based on a given large order systen. This computer 
program (OPTRED) is developed to interface with an existing 
control analysis FORTRAN prograg named OPTSYS. This thesis 
will discuss the modifications or additions to the OPTSYS 
program which are necessary to that interface. 

The OPTimal SYStems control program was originally 
developed by Hall [Ref. 1] for the study and application of 
optimal systems control theory. Later modifications tc the 
program were made by Walker [Ref. 2] and Liu [Ref. 3] of 
Stanford University. Hoden [Ref. 4] modified the OPTSYS 
program to present a user-friendly, interactive version of 
the program (OPTSYSX). Diel [Ref. 5] introduced changes 
which enabled the user to save various data for re-entry 
into the program and also to create data files which were 
necessary for the execution of atime response progran 
(OPTCALC). Purther modificaticns have been made by Laptas 
[Ref. 6] which were necessary tc create data sets for input 
to his OPTGRAPH progran. This FORTRAN program enabled the 
user to obtain Pole-Zero, Root-Locus, Nyquist, Bode and 
Nichols flots. 

It is assumed that the user has a basic understanding of 
the fundamental concepts of control theory and optinal 
systems design. The symbology conventions of Bryson 
[Ref. 7] are used in the discussion of program operation and 
system descriptions. 

An overview of the OPTSYSX progran, its capabilities, 
and the modifications to that program which are necessary 
for the interface with OPTRED is presented first. This is 
followed by the mathematical tasis upon which OPTRED is 
formulated and by a full description of the operation of 
that progran. Finally, examples of large order control 
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AND THE ACTUAL STATE #"S (IN ASCENDING ORDER) WHICH REPRESENT 
z THE REDUCED MODEL. THE ORDER OF THE REDUCED ORDER MODEL MUST 
a BE LESS THAN THE ORDER OF THE FULL SYSTEM. 


DO YOU WISH TO CONTINUE? 
TYPE "YES® OR "NO". 
DO YOU WISH RESULTS TO SCREEN OB DISK? 


NOTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARILESS OF YOUR ANSWER. 


TYPE "S" FOR SCREEN OR “D" FOR DISK. 


THE ORDER OF THE FULL SYSTEM IS: 
THE NUMBER OF CONTROLS IS: 

THE NUMBER OF OBSERVATIONS IS: 

A "D" MATRIX WILL BE INPOT. 


DO YOU WISH TO CONTINUE? 
TYPE "YES" OR "HO". 


ENTER THE DESIRED REDUCED ORDER OF THE "F" MATRIX. 


Subroutine REDUCX is now called and immediately 
reads the full system matrices from file data. The user is 
then prompted for entry of the significant states which will 
compose his reduced model. This program makes a provision 
for entry of significant states from a specifically 
formatted, user created data file. At this point in the 
program he is given that option. Once the significant 
states have been entered, the program repeats those states 
to the screen and the user is given the option to make 
changes at that time. If the significant states are satis- 
~ factory to the user, the program continues with the computa- 
tion of the reduced order nodel. At the completion of 
OPTRED, the reduced system data is created and can be exan- 
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ined for stability, controllatility and observability by 
other analysis techniques as discussed inthe following 

ee section. Prior to exiting the OPITRED program, the user is 
provided the option of re-executing OPTRED if he finds 
initial results unsatisfactory. A typical recording session 
which depicts this program flcw is presented in the next 
section. 


2. User Protection Features 


The user of OPTRED is protected from abnormal 
program termination in several ways. Initially, the progran 
presents its capabilities and user required entries. If he 
is not prepared to continue at this point, the user is 
offered the option of exiting the progran. Next, the full 
system parameters are presented and the user is once again 
able to exit if these parameters are unexpected. 

OPTRED requires that the order of the reduced systen 
ke greater than zerc but less than the order of the full 
systen. If the user inadvertantly enters a number which is 
out of this range, the program issues a warning to that 
effect and allows recovery. The computation of the reduced 
order model relies urfon the construction of a reduced state 
matrix. OPTRED requires that these states be entered in the 
order that they would normally appear in this matrix. For 
this reason, the user must enter the significant state 
numbers in ascending crder and those stute numbers must lie 
in the range from one to the rumber of states in the full 
systen. If the user inadvertantly enters these state 
humbers out of order or out of range the program frompts a 
warning and allows recovery. 

The method of reduced order computation involves the 
inversion of the state submatrix Foo as described in 
Chapter III, Equations 3.25 threugh 3.28. The nature of the 
plant matrix and the selection of reduced order states can 
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be coupled such that this matrix will be singular or naon- 
invertable. OPTRED will detect this condition and issue 
information to the user that a reduced order model cannot be 
calculated for that full system from the given desired 
reduced states. Reduced data will still be computed but the 
user is warned that this data is invalid. 

Subroutines RDINT, RDREAL and RDCHAR are responsible 
for the interactive input of data and expect an integer 
input, a real number input and a logical "YES" or "NO" input 
respectively. If the user inadvertantly makes an incompat- 
ible entry these subroutines issue warnings and allow 
another opportunity for data entry. The entry of a null 
line is included in these improper actions and the entry of 
two consecutive null strings will cause termination of the 
progran. This function allows the user a further means of 
exiting OPTRED if he so desires. 


3. Large Order Systems 


The analysis of large crder systems, particularly 
the ninety-eigth order model, presents a major problem to 
the user with limited virtual memory assets. The longitu- 
dinal control system analyzed by this thesis has a (98 X 98) 
"Pt matrix, a (2 X 98) "H" matrix and a (98 X 1) "G" matrix. 
To analyze this systen, OPTSYSX must be executed in its 
"increased dimension" form. Io compile and execute this 
version of OPTSYSX requires 2.5 Megabytes of virtual memory. 

Both OPTRED and OPTGRAPH are also dimensioned to 
enable analysis of a ninety-eigth order model. When 
OPTSYSX, OPTRED, OPTGRAPH and their peripheral data are all 
used during a typical analysis session, very careful atten- 
tion to file and virtual memory management is mandatory to 
prevent the user from exceeding his virtual machine's 
capacity. 
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C. INTERFACE WITH OTHER ANALYSTS PROGRAMS 


To sufficiently analyze a ccntrol system it is desirable 
to utilize OPTSYSX, OPTCALC, CPTGRAPH and OPTRED. It may 
also be desirable to conduct this analysis during a single, 
continuous interactive computer session. However, in the 
analysis of very large systems this concept presents serious 
difficulties in the areas of program flow, data management 
and physical memory assets. For this reason, OPTRED has not 
been developed to be fully autcmatic in its interface with 
these other analysis programs. 

The OPTRED program requires the existance of the OPTMAT 
data file. As previously discussed, OPTRED generates the 
OPTMATR data file which contains the reduced model data. To 
analyze this data using OPTSYSX, the user must now rename 
the OFTMATR data file for entry into OPTSYSX as OPTMAT data. 
Care must be exercised at this point to either erase the 
full system data file or to rename it to prevent ambiguous 
data files. Although this procedure may seen cumbersome at 
first, it is easily accomplished in the XEDIT mode and 
provides the user with the capatility of maintaining several 
system data files which are readily available for analysis. 

After OPTSYSX processes the reduced data, additional 
data files are created which enable the execution of other 
system analysis programs. OPTGRAPH and OPTCALC utilize the 
OPTMAT and OPTGROL data files, respectively, to provide the 
functions described in chapter 2. When the user is investi- 
gating several systems it is both prudent and convenient to 
rename these files after their use for the purpose of future 
7 analysis. 


De. EXAMPLES OF ORDER REDUCTION 


Two large order systems will be presented to evaluate 


reduced models generated by OPTRED. The first system is a 


Ry yaaa se Cac eter ar gia Gace ents iefayd wth ane piety eet ye teas < 
at oP eT See Oe 6 Ve Wee ae ee See ee eee pet gh ee ete eS jet Rie ater enter, ter mee to Sete eee et ea e . ee eh eee 
eh an atts, oe ea ee a ate a es wae ee at, 
Pn ne i" . ¥ we NT De er RI OO 


ve 
a er a ee a iad 
me ea a ae ee Te Same ae on Te te Me le Te Pe Me 8 


—> \ 

.. nt 
S sixteenth order model of the F100 turbofan engine. The = 
2 final example is a ninety- eighth order model of the X-29A a 
~ aircraft's longitudinal control systen. Bode analysis was ~ 
a conducted for the full and reduced models and graphical os 
es, comparisons follow at the end of the chapter. Recorded va 
Me! terminal sessions for each example will be presented to a 


fully illustrate actual program operation. 


So8 


1. Example of a Sixteenth Order System 


The F100 turbofan engine was chosen for initial 
analysis and [Ref. 9} describes this system in detail. The e 
system "F", "gn, "AH" and "Dp" matrices were obtained fron 
[Ref. 9: pp. 83-85] and the nethod of data entry to OPTSYSX 
is depicted in subroutine SETUF found in Appendix A. The eo 
selection of significant states is based upon a desired s 
control bandwidth of 1 to 10 Hertz and following modal anal- 
ysis these states were chosen as 1, 2, 5, 11 and 16. The 
following is a computer terminal system in which a fifth 4 
order model is generated. 


BEGIN RECORDING OF TERMINAL SESSION ». oe 
R; T=0.01/0.03 21:05:48 wo 
OPTRED = 
FILEDEF 05 TERM 
FILEDEF 03 DISK REDUCI DATA Al ot 
FILEDEF 06 DISK OPTRED LISTING A1 is 
FILEDEF 07 DISK STATES DATA A1 
FILEDEF 09 DISK OPTMAT DATA A1 
FILEDEF 10 DISK OPTMATR DATA A1 
GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP NONIMSL 
LOAD OPTRED ( START 
EXECUTION BEGINS... 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 23 
wp, Gy aw, AND "D" (IF INPUT) MATRICES. . 2 
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H 
THE FULL SYSTEM MATRICES MUST BE READ FROM A FILE NAMED ae 
“OPTMAT DATA" AS CREATED BY THE OPTSYSX PROGRAM. Dele 
YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) kL 


AND THE ACTUAL STATE #*'S (IN ASCENDING ORDER) WHICH REPRESENT ne 


THE RELUCED MODEL. THE ORDER OF THE REDUCED MODEL MUST BE ae 
LESS THAN THE ORDER OF THE FULL SYSTEM. pane 
be ’ - 

DO YOU WISH TO CONTINUE? eae 

as a 

TYPE "YES" OR "NO", ot 
YES 23 
DO YOU WISH RESULTS TO SCREEN OR DISK? be 

NOTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA oo 

WILL BE GENERATED REGARDLESS OF YOUR ANSWER. ies 


TYPE "S" FOR SCREEN CR "D" FOR DISK. 


s 
THE ORDER OF THE FULL SYSTEM IS: 16 r 
THE NUMBER OF CCNTROLS IS: 5 
THE NUMBER OF OBSERVATIONS IS: 7 =< 
A "D" MATRIX WILL BE INPUT. a 
Ses 
DO YOU WISH TO CONTINUE? eo 
TYPE "YES" OR "NO", hen 
YES mi 
ENTER THE DESIRED REDUCED ORDER OF THE "F" MATRIX. ete 
? aos 
; 4 
18 


*£%4*ee WARNING: REDUCED ORDER MUST BE GREATER THAN 0**#** 
AND LESS TRAN 16 


© ereele lel, 
Ship nrg mete 
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ENTER THE DESIRED REDUCED ORDER OF THE "F" MATRIX. ane 
5 = 
DO YOU WISH TO INPUT DESIRED STATES FOR YOUR REDUCED ORDER nee 
MODEL FROM A DATA FILE? r 
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DATA FILE MUST BE NAMED "STATES DATA A1" IN FIXED FORMAT. ee 
THE READ FORMAT IS "1315" PER 72 CHARACTER LINE. - 
TYPE "YES" OR "NO". eae oe 
NO Roe 
ENTER THE "N" SIGNIFICANT STATES WHICH REPRESENT THE REDUCED oe 
MODEL. ENTER STATE #'S IN ASCENCING ORDER. ke 
STATE # 1 = ee 
? = 
1 a4 
STATE # 2 = es 
? ~ 
w, 
2 oe 
STATE # 3 = : 
? : 
3 Bah 


STATE # 4 = a 


STATE # 5 = ead 
? on 
| ed 
THE REDUCED MODEL STATES ARE: oo 

1 > 3 11° «16 ee 

DO YOU WISH TO CHANGE ANY OF THE SIGNIFICANT STATES? a 

-- 

TYPE "YES" OR "NO". — 
nae hae 

a 

ENTER THE N-~TH STATE # TO BE CHANGED. as 

~*~ 

7 pe 
ENTER NEW STATE # 3 eel 

5 '. 


THE REDUCED MODEL STATES ARE: 
5 11 16 


1 2 


DO YOU WISH TO CHANGE ANY OF THE SIGNIFICANT STATES? 


NO 


TYPE "YES" 


OR "NO", 


(OUTPUT TO THE SCREEN FOLLOWS) 


THE DESIRED REDUCED ORDER IS: 5 


THE REDUCED MODEL STATES ARE: 
5 11 16 


1 2 


THE REDUCED PLANT MATRIX ("F") IS: 


-3.0509D+00 
1.2243D-01 
2.7286D-03 
&. 1314D+00 

-1.0770D+00 


THE RELUCED 


-3.5445D-02 
—-4.4094D-01 
1.7526D-02 
1.9802D+01 
-1.8909D+00 


THE RELUCED 


1.9281D-01 

6.4607D-03 

0.0 

1. 2535D-05 
-4.8870D-05 

1.6598D-05 

7.2415D-06 


re ee i a 


2-4726D+00 
-1.6271D+00 
-2.8331D-03 
-5.3854D+00 
2. 2320D+00 


-1.4080D+02 
-2.7221D+01 
~3.7646D+01 
3. 1635D+02 
7.2700D+01 


1.0750D-01 
-9.9442D-07 
0.0 
1.3512D-05 
1.4796 D-04 
-5.0355D-05 
-8.1255D-07 


eg ae Sint 


-3.58S8D+02 
4.1060D+01 
-7.4276D+00 
-1.6257D+03 
1.146€7D+03 


-9.33C1D+01 
7.8306D+00 
1.2113D-01 
4.9385D+01 

-3.0436D+01 


OUTPUT DISTRIBUTION MATRIX 


1.8753D¢02 
-9.1543D-03 

0.0 
-3.1266D-02 

1.01€0D-03 
-3.7563D-02 
-1.0493D-02 
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8.6913D+00 
4. 2898D+00 
-1. 1695D-02 
~5.3853D+01 
1. 9087D+01 


CONTROL DISTRIBUTION MATRIX ("G") IS: 


2-3597D+01 
~1.0960D+01 
~7.8465D~-02 
~7.6813D+01 
2. 8867D+01 


("H") IS: 


1. 2088D+00 
~5.6722D-05 

1.0000D+00 
~1.3900D-04 
-5. 1681D-05 
~1.7308D-04 
~4. 9442D-05 


9. 1832D-02 
3.0229D-01 
-9.4624D-03 
-2.7911D+00 
-4.8697D+01 


-1.8379D+04 
-9.9711D+03 
-1.1931D+02 
4.8330D+04 
7.0035D+03 


7.07 96D-02 
-1.9316D-05 

0.0 
-4.8476D-05 

1.4491D-05 
-5.9561D-05 
-1.6687D-05 


. 
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THE REDUCED FEEDFORWARD MATRIX 


-1.0377D-01 

-1.0929D-04 
0.0 
8.6209D-05 

-1.5078D-05 
1.0794D-04 
2.9257D-05 


3.4001D+00 
-2.281SD-01 
0.0 
-8.5310D-03 
4. 1932D-03 
5. 6356D-03 
9. 8556 D-04 


("D") IS: 
1.0378D+01 1.7019D+00 
3.1109D-01 4.8490D-03 
0.0 0.0 

-4,4940D-03 -9.5892D-05 

-9.3797D-04 -3.9608D-03 

-4.0947D-04 -7.9878D-04 
3.2771D-04 -9.4642D-05 


-4.1189D+03 
-1.31705D+01 
0.0 
~4.4513D-01 
9.6895D-01 
-3. 1434D-01 
-8.3787D-02 


ANALYSIS COMPLETE... YOUR REDUCEL SYSTEM DATA HAS BEEN SAVED 
IN A FILE NAMED “OPTMATR DATA". 


NO 


IF YES, 


DO YOU DESIRE ANOTHER RUN? 


THE RESULTS OF THIS RUN WILL 


REPLACE YOUR LAST OPIMATR DATA. 


TYPE "YES" 


CR "NO", 


eeeveeeeeOPTRED IS NOW TERMINATED... -.. 0 


R; T=0.38/0.77 21:06:56 


RECORD OFF 


END RECORDING OF TERMINAL SESSION 


This is a typical computer terminal session in which 


tke user requested output to 
some of the user protection 


during this session. 


Following this analysis, the 


renamed OPTMAT data 
reduced model data. 

able characteristics, including stability, 
to the full system eigenvalues. 
factory and further 


and OPTSYSX 
The eigensystem was analyzed for desir- 


comparison 
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the screen. 


Note also, 


that 


features of OPTRED are involked 


OPTMATR data file was 


was executed 


and 


using the 


and was compared 
This analysis proved satis- 
of the full 


reduced 
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models was accomplished by obtaining Bode diagrams for the 
system ofen loop transfer functions. Figures 4.1 and 4.2 
illustrate the full system and reduced system Bode plots 
respectively. The similarity, especially within the control 
bandwidth, is readily apparent. The frequency responses 
shown are for system input 1 and output 1. Analysis of 
other input/output combinations revealed the same favorable 
comparisons. 


2. Example of a Ninety-eighth Order Systea 


The X-29A control system includes a ninety-eiyhth 
order dynamics matrix. The measurement scaling matrix is 
constructed such that there are two measured outputs. These 
outputs represent the system with and without a notch 
filter, respectively. Full system matrix data as well as 
eigensysten data was obtained for this model fron 
NASA-Edwards. Full system data was first entered into 
OPTSYSX using subroutine SETUP and the overall system was 
analyzed. Sixty significant states were identified and are 
shown in the following recorded terminal session. This 
example also illustrates program execution where no "“p" 
matrix is input and the reduced system states are entered by 
means of file data. During this session the user also 
chooses to have the general reduced model data printed toa 
listing file. 


BEGIN RECORDING OF TERMINAL SESSION 
R; T=0.01/0.02 12:23:56 
OPTRED 
FILEDEF 05 TERM 
FILEDEF 06 DISK OPTRED LISTING Al 
FILEDEF 07 DISK STATES DATA At 
FILEDEF 09 DISK OPTMAT DATA Al 
FILEDEF 10 DISK OPTMATR DATA A1 
GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP NONIMSL 
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LOAD OPTRED ( START 
EXECUTION BEGINS... 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
NPA GH UHM, AND "D" (IF INPUT) MATRICES. 


THE FULL SYSTEM MATRICES MUST BE READ FROM A FILE NAMED 
“OPTMAT DATA" AS CREATED BY THE OPTSYSX PROGRAM. 


YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) 
AND THE ACTUAL STATE #*S (IN ASCENDING ORDER) WHICH REPRESENT 
THE RELTUCED MODEL. THE ORDER OF THE REDUCED MODEL MUST Be 

LESS THAN THE ORDER OF THE FULL SYSTEM. 


DO YOU WISH TO CCNTINUE? 


TYPE "YES" OR "NO", 


YES 
DO YOU WISH RESULTS TO SCREEN OR DISK? 
NCTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARDLESS OF YOUR ANSWER. 
TYPE "S" FOR SCREEN OK “"D" FOR DISK. 
D 
THE ORDER OF THE FULL SYSTEM IS: 98 
THE NUMBER OF CONTROLS IS: 1 
THE NUMBER OF OBSERVATIONS IS: 2 
A "D" MATRIX WILL NOT BE INPUT. 
DO YOU WISH TO CONTINUE? 
TYPE "YES" OR "NO". 
YES 
ENTER THE DESIRED REDUCED CRDER OF THE "F" MATRIX. 
? 
60 


DO YOU WISH TO INPUT DESIRED STATES FOR YOUR REDUCED ORDER 
MODEL FROM A DATA FILE? 
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DATA FILE MUST BE NAMED "STATES DATA A1" IN FIXED FORMAT. nee 
THE READ FORMAT IS "13I5" PER 72 CHARACTER LINE. 
TYPE "YES" OR "NO". 
YES 


THE KEDUCED MODEL STATES ARE: 

1 #2 3 4&8 § 6 42 #13 4% 1 43 44 45 i 

4647 #48 49 5SO 51 52 53 54 55 62 63 64 

65 66 67 68 69 7 71 #72 #73 «74 #«+IS 7% 77 

73 #79 #80 81 82 83 84% 85 86 87 88 89 90 > 
91 92 93 94% 95 96 97° 98 


DO YOU WISH TO CHANGE ANY OF THE SIGNIFICANT STATES? 


TYPE "YES" OR "NO", > 
NO 


ANALYSIS COMPLETE... YOUR REDUCED SYSTEM DATA HAS BEEN SAVED 
IN A FILE NAMED "OFTMATR DATA". 


DO YOU DESIRE ANCTHER RUN? 


IF YES, THE RESULTS OF THIS RUN WILL 
REPLACE YOUR LAST CPIMATR DATA. 
TYPE "YES" OF "NO", he 
NO : 


eeeceee-OPTRED IS NOW TERMINATED........ 4 


Rs; T=17.70/18.51 12:25:07 
RECORD OFF 
END RECORDING OF TERMINAL SESSICN 


The eigensystems for both systems were analyzed for 


stability and other desireable characteristics. Bode piots 


were obtained for both input/output combinations and these 
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gtaphs for the full and reduced systems are included as 
Figures 4.3 through 4.10. The input/output combination 
(1/1) shows the response of the system without the notch 
filter while the (1/2) input/output combination depicts the 
frejguency cesponse of the system after the addition of the 
notch filter. These figures clearly indicate excellent 
agreement between the full and reduced order systems in the 


frequency domain. 
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thesis. The addition of this yrogram to the analysis tech- oS 
nique will provide a further means of verification of the = 
accuracy and fidelity of the reduced order model. 


2. Interface wit 


The Naval Postgraduate School's Aeronautical 
Engineering Department has developed a collection of anal- 
ysis programs including OPTSYSX, OPTGRAPH and OPTCALC. 
These programs exist ina single "module" which enables the 
user to interactiveiy execute these programs during a single 
computer terminal session. The feasibility of adding OPTRED 
to this rodule should be investigated. If necessary, these 
programs would require re-dimensioning and could create a 
very large memory size reguirement. Additions and modifica- 
tions would also be necessary in the OPTSYSX EXEC to enable 
the execution of OPTRED and then transfer execution to ae 
OPTSYSX for the analysis of the reduced model. This 
controlling EXEC must also manage the data files upon which 
the individual programs depend. 


“few ee 
a ea ee 


3. Alternate Order Reduction Techniques eS 


Pinally, other methods of order reduction should be a 
investigated. Comparisons of reduced order models which are : 
generated by different techniques would provide the controls ie 


designer with additional means cf reduced model verification 
and analysis. 
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A. CONCLUSIONS 


The computational capabilities of OPTRED were tested by 
analyzing two significant large order systems. The deterni- 
nation of the validity of both full order and reduced order 
results was greatly simplified by the eigensystem data for 
the X-29A provided by NASA-Edwards. As discussed by Laptas oo 
[Ref. 6], OPTSYSX failed to produce valid open-loop transfer 
function zeroes when analyzing very large systems. However, 
through modal analysis and with the aid of the NASA data, 
the extraneous zeroes were easily identified and the open L 
loop data created by OPTSYSX was corrected. 

The results of OPTSYSX and OPTGRAPH indicate a very i 
strong similarity between the full system and the reduced ce 
system generated by OPTRED both in stability and in the 


frequency domain. The combined analytic capabilities of oe 
OPTSYSX, OPTGRAPH and OPTRED provide a powerful tool for the pe 
controls designer for use in the analysis, development and ad 
the implementation of complex large order systems. — 


B. RECOMMENDATIONS 


Based upon the results of this thesis the following 
areas are recommended for further study and research. —- 
1. Comparison of System Time Response 


Although the reduced model must closely resemble the 
full system in the frequency domain, it should also compare 
favorably in the time domain. OPTCALC, the time response 
program, requires several modifications to enable its use 
for analysis of systems as large as those presented in this 
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